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ABSTRACT 



We investigate the effect of dust on spiral galaxies by measuring the inclination- 
dependence of optical colours for 24,276 well-resolved SDSS galaxies visually classified 
via the Galaxy Zoo project. We find clear trends of reddening with inclination which 
imply a total extinction from face-on to edge-on of 0.7, 0.6, 0.5 and 0.4 magnitudes 
for the ugri passbands (estimating 0.3 magnitudes of extinction in z-band). We split 
the sample into "bulgy" (early-type) and "disky" (late-type) spirals using the SDSS 
f racdeV (or foev) parameter and show that the average face-on colour of "bulgy" 
spirals is redder than the average edge-on colour of "disky" spirals. This shows that 
the observed optical colour of a spiral galaxy is determined almost equally by the spiral 
type (via the bulge-disk ratio and stellar populations), and reddening due to dust. We 
find that both luminosity and spiral type affect the total amount of extinction, with 
disky spirals at M r ~ —21.5 mags having the most reddening - more than twice as 
much as both the lowest luminosity and most massive, bulge-dominated spirals. An 
increase in dust content is well known for more luminous galaxies, but the decrease of 
the trend for the most luminous has not been observed before and may be related to 
■ their lower levels of recent star formation. We compare our results with the latest dust 

attentuation models of Tuffs et al. We find that the model reproduces the observed 
trends reasonably well but overpredicts the amount of u-band attenuation in edge-on 
galaxies. This could be an inadequacy in the Milky Way extinction law (when applied 
to external galaxies), but more likely indicates the need for a wider range of dust- 
star geometries. We end by discussing the effects of dust on large galaxy surveys and 
emphasize that these effects will become important as we push to higher precision 
measurements of galaxy properties and their clustering. 

Key words: galaxies: spiral - galaxies: fundamental parameters - galaxies: photom- 
etry - ISM: dust, extinction - surveys 
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1 INTRODUCTION 

The clear view we enjoy of the extragalactic sky towards the 
Galactic poles led to an early assumption that most disks of 
spiral galaxies are largely transp arent. Early work supported 
this i dea (e.g. [H olmbcrg] ll958h . and it was not challenged 
until iDisnev et all (|lQ89h pointed out that the FIR emis- 
sion observed by the IRAS satellite could not be explained 
without absorption and re-emission of optical light by sig- 
nificant amounts of dust. Since then several authors have 
used the inclination dependence of dust extinction to judge 
the role of dust on the observed properties of spiral galaxies 
fee. ICunowHl992l: iGiovanelli et all Il994l : iTullv et all 1 19981 : 
iMasters et al.ll2003h . and it remains unclear even today if 
most spiral galaxies are predominantly optically thin (low 
opacity) or optically thick (high opacity). Studies of over- 
lapping galaxies have shown that the relative geometry of 
the stars and dust plays an important role, which can vary 
significantly from galaxy to galaxy. They also show that in 
a given galaxy there may be both optically thick and thin 
regions which may or may not correlate w ith patterns in the 
stellar density (e.g. iHolwerda et al.ll2009T ). 

In recent years, several authors have revisited this prob- 
lem using the sheer size and quality of la rge galaxy survey s 
like the Sloan Digital Sky Survey (SDSS; lYork et alj|200cj ). 
For example, Alam & Ryden (2002) first pointed out that 
the red population of SDSS galaxies (selected via a u — r 
color cut) is contami n ated by dust-reddened edge-on spi- 
rals, while Sha o et all (120071 ^) studied the dependence of the 
luminosity function of 61,506 SDSS spiral galaxies (selected 
using /DevEb: 0.5) on inclination, finding that dust extinc- 
tion caused about 0.5 magnitudes of dimming in z-band and 
1.2 mags in u-band, cons istent with what would be e xpected 
for optically thick disks. lUnterborn &: Rvdenl (|2008l ) used a 
more stringent cut of /d c v < 0.1 (arguing that fuev < 0.5 
will result in many early- type interlopers) to select 36,162 
late-type spirals. They study trends of u — r colour and r- 
band magnitude in a volume limited subset and find ~ 1.3 
mags of dimming from face-on to edge-on in r-band. 

These inclination effects on the global properties of 
galaxies have been discussed in most detail by Driver et 
al. (2007) usi ng data from th e Millenium Galaxy Cata- 
logue ( MGC | Liske et all I2003T ) wit h bulge-disk decompo- 
sitions |Allen et al. 1120061 ) as well as lMaller et alj l|2009h us- 
ing SDSS data. Both of these studies highlight that most 
measured galactic distributions and relationships, especially 
from the SDSS, are biased by dust effects which can be up 
to 2-3 magnitudes in shorter wavelength bands for the most 
inclined galaxies. Even the original target selection for the 
SDSS Main Galaxies (Strauss et al. 2002), made in r-band 
could be affected by these issues and this likely leads to sub- 
tle incompletenesses in the studies of the large-scale cluster- 
ing of galaxies (since elliptical galaxies and spiral galaxies 
have different clustering properties). 

We have revisited this issue using a new sample of spiral 



1 The SDSS measures the fracdeV or /d c v parameter, which 
describes the fraction of the light fit by a de Vaucouleurs profile 
verses an exponential profile. A pure de Vaucouleurs elliptical 
should have /o e v = 1> and a pure exponential disk spiral will 
have / DcV = 0. 



galaxies selected from the Galaxy Zoo project^. All previous 
SDSS studies of the inclination effects have used some mea- 
sured proxy (color, model fits, concentration, Sersic index) 
to select their spiral or disk galaxies. As discussed above, 
there is much debate over the best parameter values to use in 
selecting such galaxies as well as significant scatter between 
these different proxies (see Appendix A). In contrast, Galaxy 
Zoo provides robust visual classifications for over 10 5 objects 
thanks to the partic ipation of more than 160,000 volunteers 
(|Lintott et alj 12008), and in recent Galaxy Zoo papers we 
have demonstrated the complexity of relating galaxy colors 
(and other mor phological parameters ) to these visual clas - 
sifications, e.g. , 
IMasters et al. I 



Bamford et all feOOSMSkibba et al l i|2009f) ; 
201Ch , and lSchawinski et al. I (|2009l 'l who dis- 



cuss the interesting sub-populations of "red spirals" and 
"blue ellipticals" respectively, found in the Galaxy Zoo data. 

In this paper, we study the trends for A — z (ie. op- 
tical colours relative to the z-band) reddening as a func- 
tion of axial ratio for a sample of well-resolved Galaxy 
Zoo spiral galaxies. While the total e xtinction in z-band 
is n ot expected to be z ero (for example IMasters et al.l [20031 
and [5 river et al. I l2008l both show non-zero K-band extinc- 
tion), we avoid comb ining SDSS and 2MASS (or UKIDSS, 
lLawrence et al.l|2007l ) data due to worries about systematic 
differences in the photometric apertures between the SDSS 
and these near-Infrared (NIR) surveys, as a function of in- 
clination. In the future, such studies would be improved by 
including the NIR data and thus giving a "zero extinction" 
measurement. 

In Section 2, we describe the data from Galaxy Zoo 
and SDSS, including galaxy photometry, and axial ratios, 
and discussed biases introduced by our sample selection. In 
Section 3, we show the observed change in color as a function 
of axial ratio, and then in Section 4, compare these t rends 
to the dust attenuation model from iTuffs et al. 1 1 |2004 ). We 
discuss the implications of the work and conclude in Section 
5. 



2 DATA AND SAMPLE SELECTION 
2.1 Galaxy Zoo Classifications 

The Galaxy Zoo (GZ) project uses an internet tool to al- 
low volunteers from the general public to visually classify 
galaxies observed by SDSS (see iLintott et al.l 120081 for de- 
tails of the sample selection and initial results). In particu- 
lar, Galaxy Zoo participants were asked to say if a galaxy 
was elliptical, spiral, "don't know" or a merger. The spi- 
ral classification was then divided into either clockwise or 
anti-clockwise classes (based on the apparent direction of 
the spiral arms), or "edge-on/don't know". In total, each 
SDSS galaxy received an average of 38 separate classifica- 
tions, with most of the galaxies having at least 20 inde- 
pendent classifications. Based on these classifications, each 
object is then assigned a likelihood of being either a spi- 
ral or an elliptical (or merger or "don't know"). Results 
using these Galaxy Zoo class ifications have been presented 



in a series of recent papers 1 Land et al.l 



120091 ; iBamford et all 120091 ; ISkibba et al 



20081; ISlosar et all 



20091 ; iDarg et all 
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2009allbT: ISchawinski et all 120091 : ICardamone' et al. 2009; 



Masters et al~ll201Ch 



The Galaxy Zoo project was initiated by the need for 
reliable visual morphologies of SDSS galaxies - a sample an 
order of magnitude larger than any which had previous vi- 
sual classifications. Discussions of why automated methods 
for classification were deemed insufficient for many scien- 
tiflc purposes can be f ound in the introduc tions of both 
iLintott et ail <|2008l ) and lBamford et all <|2009h . We show, in 
Appendix A, that visual classification is important to avoid 
the high levels of incompleteness and contamination present 
in morphological samples selected using simple structural 
measurements. We define a sample of spiral galaxies us- 
ing the Galaxy Zoo classification probabilities, p sp irai > 0.8. 
We include a redshift cut of O.Of < z < 0.09 to reduce 
the redshift cla ssification bias discussed in Appendix A of 
Bam ford et al.l (120091'). and we have als o applied the correc- 
tions outlined in lBamford et all {2009). These cuts provide 
a sample of 79,935 visually classified spiral galaxies, which 
to M r — —20.5 is a volume limited sample for face-on spirals 



2.2 SDSS Photometry 

Our photometric quantities are taken from the S DSS Data 
Release 6 fDR6. lAdelman-McCarthv et all l2008h . We use 
Petrosian magnitides for the total flux, and model magni- 
tudes for colours. These are corre cted for Galactic e xtinc- 
tion using the DIRBE dust maps (|Schlegel et alJll998h and 
have a small k-corr e ction applied using KCO RRECTv4_f _4 
(|Blanton et alj|2003t iBlanton fc Roweid l2007f ). Throughout 
this paper we use a standard cosmology (ft&f = 0.3, fl\ = 
0.7) with H = 70km s" 1 Mpc" 1 . 

As a measure of the size of galaxies, we use the radius 
enclosing 90% of the Petrosian flux (in r-band), rgo. We also 
use the concentration index, c = rgo/rso (the ratio between 
the radius enclosing 90% of the Petrosian flux and 50% of 
the Petrosian flux). We remind the reader that a Petrosian 
radius is defined to be the radius where the local surface 
brightness (ie. the mean surface brightness in a small an- 
nulus) is equal to a constant f raction of the mean surface 
brightness within that radius (jStrauss et al J [2002? ) . fn the 
absence of dust, or in a completely opaque disk, the Pet- 
rosia n radius is independent of inclination fsee lStrauss et ail 
2002). Dust will cause the light profile of a galax y to become 
shallower (as described in lGiovanelli et ailll994l ') causing the 
Petrosian radius , and rgo, to move outwards. For example 
iMollenhoff et al] (|2006l ) show that in B-band sizes can in- 
crease by f 0-40% from face-on to edge-on in the presence of 
dust. 

We also use the SDSS structural parameter f racdeV (or 
/uev)- ln SDSS, all galaxy light curves are fit with both a de 
Vaucouleurs and an exponential profile. Model magnitudes 
are contructed from an optimal combination of these two 
fits (the linear combination which best fits the data), and 
fDeV describes the fraction of the light which is fit by the 
de Vaucouleurs profile. 



inclination will vary from galaxy to galaxy. It will depend 
on the intrinsic axial ratio, q, of the galaxy (i.e. that which 
would be measured for i = 90°) and on any variation from 
a/b — 1 which would be observed if the galaxy were face- 
on. A reasonab le intrinsic axial ratio for spirals varies from 
q = 0.1 to 0.2 jllnterborn fc Rvdenll2008l measure q ~ 0.22), 
while the ellipticity of the face-on disk is usually assumed to 
be small (it has been meas ured at e ~ 0.08, or b/a — 0.92 by 
lUnterborn fc Rvdenlliooj ) so is neglected. Inclination can 
therefore be estimated using 



{b/af-q 2 
I - q 2 



(1) 



We use the p-band isophotal a/b reported by SDSS 
which are fit to the 25th mag/square arcsecond isophote 
(t"2s)- We pick this axial ratio rather than that calculated 
from the flux weighted moments of the galaxy as it is likely 
to be a better estimate of the maximum ellipticity and there- 
fore of the true inclination of the galaxy (since the moments 
a/b is flux weighted it is likely to reflect the shape in the 
brightest parts of the galaxy) . We pick g-band as a compro- 
mise between getting the best S/N (which would lead us to 
pick the default r-band) and the fact that bluer bands better 
trace the shape of the disks of spirals at these redshifts. 

Obviously, both the bulge and disk components of a spi- 
ral affect it's observed axial ratio. The best way to measure 
the inclination is through the axial ratio of the disk in a 
bulge-disk decomposition, but in the absence of that data a 
value of q in Eqn (I) which increases with B/T provides the 
best possible estimate. 

As discussed on the SDSS Algorithms Pag^f]. the 
isophotal axial ratio is not corrected for seeing. We see the 
impact of this in that the axial ratios of the smallest GZ spi- 
rals are on average rounder than the larger spirals. Figure [T] 
shows rgo against log(a/fo) measured in g-band. Horizontal 
lines are plotted at rgo = 4" (the apparent minimum of this 
value for GZ spirals at O.Of < z < 0.09) and r go = 10", 
where the effects of seeing appear to disappear. The curved 
line shows a simple model in which the minimum possible 
semi-minor axis is b — 1 .3" (an estimate of the mean see- 
ing in SDSS), and the major axis is rgo. This appears to 
match the observed trend well. It is possible to correct the 
ob served axial ra t ios fo r the effects of seeing as was done 
by iMasters et ail {2003) for a sample of 15,224 spirals in 
2MASS. However, given the size of our sample, we choose to 
deal with this issue by simply selecting only large GZ spirals 
where seeing is not an issue, i.e., rgo > 10". In the rest of 
this paper, we refer to this selection as our "well resolved" 
GZ spirals (O.Of < z < 0.09) which totals 24,279 galaxies. 
This selection still admits galaxies over the full redshift and 
luminosity range above. A Milky Way size galaxy might just 
make it into the selection at the high redshift limit where f 0" 
corresponds to a physical size of about 17 kpc (it's roughly 
2 kpc at the lower redshift limit). However as discussed be- 
low in Section 2.5, this strict cut on size will complicate the 
biases in our sample. 



2.3 Inclinations 

We use the observed axial ratio (a/b, from SDSS) as a proxy 

for inclination. The exact correspondence between a/b and 3 www.sdss.org/dr7/algorithms/ 
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Figure 1. The g-band axial ratio versus the radius enclosing 
90% of the r-band light for GZ spirals (p spiral > 0.8, 0.01 < z < 
0.09). This shows the effects of seeing which rounds the isophotes 
of small galaxies. Horizontal lines are plotted at rgo = 4" (the 
apparent minimum of this value for GZ spirals at 0.01 < z < 0.09) 
and rgo = 10", where the effects of seeing appear to disappear 
almost completely. The curved line shows a simple model in which 
the minimum possible semi-minor axis is b = 1.3" (an estimate 
of the 25th percentile best seeing in SDSS - the median is about 
1.43"). 

2.4 Spiral Type or Bulge-Disk Ratio 

As has been known since iHubbld l|l926l ) not all spiral 
galaxies are the same. Both their physical properties and 
appearance vary significantly from the "early type" spi- 
rals (Sas with large bulges) to the "late type" spirals 
(Sc/Sds with little / no bu lge). For a review of this topic see 
iRoberts fe Havnesl (|l994l ). Since the star formation histories 
of different types of spirals differ (with later types having 
more recent star formation), so do their intrinsic colours. 
The expected intrinsic axial ratio (how thick they appear 
when viewed edge-on) is also affected by the presence of a 
bulge, such that earlier type spirals with large bulges are 
wider (have smaller a/b) when viewed completely edge-on 
than their late type spirals counterparts. 

Therefore, to fully understand the reddening of spirals 
with their observed axial ratio we need a way to split galax- 
ies of differing bulge-disk ratios and intrinsic axial ratios. 
This has been done using bulge-disk decompositions to study 
the B-band attenua tion of 10095 galax ies in the Millennium 
Galaxy Catalogue (|Driver et al.ll2007h . however SDSS does 
not provide a bulge-disk decomposition as a standard pa- 
rameter. We therefore, attempt to use other structural pa- 
rameters from SDSS to broadly divide the sample into those 
spirals with a pure disks (i.e. late type spirals) and those 
with large bulges (i.e. early type spirals). 

2.4-1 Finding Bulges with Light Profile Information 

In order to split the sample by "bulginess" we use struc- 
tural parameters provided by SDSS, namely concentration, 
c, and /oev- A classic elliptical galaxy is expected to have 
values of c ~ 5.5, and /oev = 1 while a pure exponen- 
tial d isk will have c ~ 2.3 and fDeV = (|Strateva et al.l 
2001). Based on a sample of ~ 300 visually classified galax- 



ies, IStrateva et al.l |200ll ) recommend that c = 2.6 be used 
to divide the population. Furthermore, they argued that a 
galaxy with fDeV > 0.5 is likely to be an elliptical, SO or 
Sa galaxy, while those with /oev < 0.5 are late type spirals 
(Sb or Sc) or irregular. (We explore these cuts with Galaxy 
Zoo classifications in Appendix A.) 

In Figure [2] we plot both concentration and fDeV as 
a function of axial ratio and (u — z) colour for our sample 
of "well resolved" GZ spirals. It is clear that concentration 
is not a clean method for spliting our spirals as there is a 
strong trend such that more inclined spirals are more likely 
to have high concentrations. The best linear fit to this trend 
is c = 2.169(3) + 0.632(9) log(o/6) 0. Even the most inclined 
spirals do not reach the c = 5.5 expected for a classic el- 
liptical galaxy, however they do reach the lowest concentra- 
tions of the observed elliptical range which as discussed in 
IStrateva et al.l l|200ll) has a large scatter. The mean value 
of the concentration of Galaxy Zoo spirals reaches c = 2.6 
at log(a/b) ~ 0.7 or i ~ 80° for late type spirals (assumes 
q — 0.1), so 50% of these edge-on spirals will fall into the 
"early type" subset if such a divider is used (also see Ap- 
pendix A). 

On the other hand fDeV appears to be a good candi- 
date to split the spirals into t hose with large bulge s, and 
those with no bulge. Ironically, IStrateva et alj f2001) warn 
about the use of foev to split galaxies by type in large, 
bright galaxies as the model is dominated by light from the 
central regions, but herein, this is exactly the signal we are 
looking for, i.e., spirals with large central bulges. We see 
in Figure [2] that "blue" spirals are found to be only those 
with small fneVi while "red" spirals span the full range of 
fDeV which we interpret as a mix of intrinsically red face-on 
spirals with large bulges and dust reddened edge-on spirals 
with no/small bulge. Supporting this interpretation is the 
fact that there are no galaxies with both large values of 
log (a/b) (ie. very thin disks viewed edge-on) and large val- 
ues of fDeV, (i-e., large bulges). 

As further evidence that the fDeV parameter is able 
to divide the visually classified GZ spirals by Hubble type 
(or bulge size), we show in Figure [3] the colours of face- 
on (log(a/6) < 0.05) GZ spirals plotted against their /zjev- 
As expected for early-type spirals (with large bulges), the 
objects with the largest values of fneV always have redder 
(mean) face-on colours than those with small values of fDeV, 
(i.e., spirals with smaller bulges). A t first glance this might 
seem to contradict the findings of iDrorv fc Fisheij (|2007l ) 
who show in a sample of 39 objects that disk galaxies with 
classical bulges are globally red (irrespective of the bulge 
size), while only in galaxies with pseudo bulges is there a 
trend of global colour with B/T (at low values of B/T - 
since pseudobulges are not found at high B/T). This con- 
tradiction cannot be clearly tested without identifying the 
classical and pseudo- bulges in our sample face-on spirals, 
however looking at Figure [3] in more detail we argue that it 
may still show the expected trend. There is a larger spread 
of colours of our face-on spirals with small values of fz>eV 
than is seen in those with large values of fDeV - consistent 



4 We are using a condensed notation for quoting errors on fitted 
parameters where 2.169(3) is equivalent to 2.169 ± 0.003. This is 
used throughout the paper. 
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Figure 2. Concentration and foeV f° r our sample of well resolved GZ spirals as a function of observed axial ratio and (u — z) colour. We 
plot the best fit to the relation between concentration and axial ratio, and also a fit to our estimate of the relation between the intrinsic 
axial ratio (ie. the maximum value of a/b possible) with foeV- 



with the idea of a mix of redder galaxies with classical bulges 
and bluer galaxies with pseudo-bulges at low fDeV and only 
redder galaxies with classical bulges at high values of fDeV- 

We note that our data shows that using a strict cut on 
fpev (as done bv lShao et al.|[2007l and lUnterborn fe Rvdenl 
2008) will miss a large fraction of the spiral population. It 
also shows that an early-type sample selected using a mini- 
mum value of fDeV will have significant contamination from 
visually classified spirals with large values of fDeV- We dis- 
cuss this further in Appendix A. 

In Figure SI we show the match between our well- 
resolve d GZ spirals and the Millennium Galaxy Catalog 
(MGC, iLiske et al.ll2003h , finding an overlap of 109 galaxies 
(the MGC is much deeper than SDSS but over a smaller 
area). We use this match to find a linear relation between 
the /oev parameter from SDSS and the bulge fraction as 
measured by the MGC. This relation is shown in Figure [4] 
and is given by f De v = 0.11(3) + 1.7(2)(S/T)mgc 




Figure 4. SDSS foeV parameter against the MGC bulge-fraction 
for the 109 well-resolved GZ spirals found in the MGC. The solid 
line shows the best fit value of f DeV = 0.11(3) + 1.7(2)B/T. 



Finally, Figure [5] shows example images of face-on GZ 
spirals of similar angular size, ordered by foeV- 
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Figure 3. The colour of face-on GZ spirals (log(a/fe) < 0.05) plotted against foeV- The solid line in each panel is a linear fit to the 
data, and the dots show the medians in 17 bins of 100 galaxies each. 



2.4.2 Intrinsic Axial Ratios 

In order to obtain the inclination of a galaxy from its 
observed axial ratio, an estimate of its intrinsic axial ra- 
tio (ie. the axial ratio which would be observed if it were 
completely edge-on) is needed. This is observed to vary 
along the Hubble sequence, with earlier type spirals with 
large bulges app earing thicker when viewed edge on than 
later type spirals. lHavnes Sz Giovanellil |l984]) tabulated the 
measured values of q for different Hubble types (based on 
work from an unpublished prep rint by Lewis 1980), while 
ISimien fc de Vaucouleursl |l98a) give the expected bulge- 
total luminosity ratios for different Hubble Types. We pro- 
vide in Table [TJ a compilation of these two results in which 
we also include the approximate value of fn e v we find for 
these values of B/D in our well resolved GZ spirals. 

Here, we use foeV to predict the intrinsic axial ratio 
of spiral galaxies. We plot in the bottom left-hand panel 
of Figure [2] an estimate of the maximum axial ratio as a 
function of fDeV for which we use log(a/&) max = (log(a/fe)) + 
2o"io g (a/b) ■ We will use a fit to q = (fe/a) max from these values 
as an estimate of the trend of the intrinsic axial ratio of a 
galaxy with foeV, which we find to be q = 0.12 + 0.10/deV, 
consistent with the range of value of q observed in edge-on 
galaxies in Table [TJ 



Table 1. Bulge-Disk Ratio and Intrinsic Axial Ratio as a Func- 
tion of Hubble Type. Also included is an estimate of the mean 
value of fDeV f° r each type. 



Type 


T 


B/T 


B/D 


q 


fDeV 


SO- 


-3 


0.6 


1.5 


0.23 


1.0 


SO 


-2 


0.57 


1.3 


0.23 


1.0 


S0+ 


-1 


0.53 


1.1 


0.23 


1.0 


SO/a 





0.48 


0.92 


0.23 


0.9 


Sa 


1 


0.41 


0.69 


0.23 


0.8 


Sab 


2 


0.32 


0.47 


0.23 


0.6 


Sb 


3 


0.24 


0.32 


0.23 


0.5 


Sbc 


•1 


0.16 


0.19 


0.2 


0.4 


Sc 


5 


0.094 


0.10 


0.175 


0.3 


Scd 


6 


0.049 


0.05 


0.14 


0.2 


Sd 


7 


0.022 


0.022 


0.103 


0.1 



2.5 Sample Bias 

Studies of the trends of observed quantities with inclina- 
tion are hampered by sample bias. Most observed quantities 
depend in some way on inclination, therefore cuts on ob- 
served quantities produce samples biased by inclination. The 
intercorrelations between properties of galaxies (for exam- 
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Figure 5. Examples of face-on GZ spiral galaxies ordered by foeV 



pie size-magnitude and colour-magnitude relations) also add 
complications. An unbiased sample for the study of inclina- 
tion trends is probably not possible, but an appreciation of 
the effect of any bias in the sample on the results will aid in 
the interpretation. 

While the majority of face-on SOs will not make it into 
the GZ spiral cut, even expert classifiers are more likely to 
classify SOs as spirals if they are observed edge-on. Therefore 
we expect the edge-on GZ spirals may have a larger propor- 
tion of SO galaxies than the face-on. A comp arison of ex- 
pert classifications from Fukugita ct al. I (I2007T) and Galaxy 
Zoo classifications (|Lintott et alJl2008l : iBamford et al.ll2009 ) 
shows tha t SOs contribute less t han 3% to the GZ spirals, 
however if iFukugita et all (|2007h also tend to classify edge- 
on SOs as spirals this comparison is not useful. The pos- 
sible "contamination" of SOs into the edge-on GZ spirals 
should be remembered when interpreting our results - since 
SOs are usually brighter and redder than spiral galaxies this 
could cause attenuation to be underestimated while redden- 
ing measures could be overestimated. 

The sample selected here, based on the SDSS Main 
Galaxy Sample has an implicit magnitude limit in the r- 
band of m r = 17.77. Our redshift limit of 0.01 < z < 0.09 
therefore makes this a volume limited sample for face-on spi- 
rals to M r = —20.5, while dimmer spirals are only observed 
in the near parts of the sample. Because of the dimming 



effects of dust, any magnitude limited sample will miss in- 
trinsically brighter objects at large inclinations. 

We also apply a size cut to our sample. (As discussed 
above, because of the effects of seeing, we cut at an ob- 
served Petrosian 90% flux semi-major axis (in r-band) of 
?~9o = 10".) The presence of dust in a galaxy disk shallows 
the light profile and moves 7-90 outwards as a galaxy becomes 
more inclined. This size cut then allows intrinsically smaller 
objects into the sample at larger inclinations. There is an ob- 
served correlation between size and magnitude (most small 
galaxies are also dimmer) however there is also scatter in 
that relation, so the net effect here is to have a sample of 
inclined galaxies which is intrinsically brighter and smaller 
(i.e. more compact) than the comparison face-on galaxies. 

We show in Figure [S] the observed size-magnitude rela- 
tion for our sample of well resolved GZ spirals. As is well 
know, brighter galaxies are physically larger. As expected, 
the most inclined galaxies in our sample have observed sizes 
which are larger than the more face-on at a given observed 
(k-corrected and Galactic extinction corrected) magnitude. 
We expect that intrinsically these two populations should 
follow the same relation, so this illustrates the effect of in- 
clination on our sample selection - it dims the observed 
magnitudes and increases the radii outwards at the same 
time. To make these two relations align, the magnitudes of 
the edge-on spirals must be shifted by ~0. 5-1.0 magnitude 
(ignoring any change in size), quite consistent with the re- 
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Figure 6. Observed physical size versus absolute r-band mag- 
nitude (corrected for Galactic dust and k-corrections, but not 
internal extinction) for our sample of well resolved GZ spirals. 
The most face-on quartile are hi-lighted in blue, while the most 
edge-on quartile are hi-lighted in red. 
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Figure 8. Estimate of the cos i distribution of the full sample. 
Poisson counting errors are shown, as is the expected flat distri- 
bution (dashed line). The estimated median error on the value of 
cos i is shown by the horizontal bar. 
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Figure 7. Absolute r-band magnitude (upper panel) and phys- 
ical size (lower panel) of the sample of well resolved GZ spirals 
discussed in this paper. Only the most face-on (blue points) and 
most inclined (red points) quartiles are shown to illustrate the 
differences in the samples at the two extremes of inclination. 

suits on total r-band extinction found here (see Section 3) 
and other studies (e.g. lUnterborn Sz Rvdenll200ot ). Figure [7] 
shows the result of the two cuts (SDSS Main Galaxy Sam- 
ple magnitude limit, and our applied radius limit) on the 
absolute magnitude and linear sizes of the sample as a func- 
tion of redshift. Only the most face-on and edge-on quar- 
tiles of the sample are shown to illustrate the differences. At 
all redshifts the observed magnitudes of the inclined galax- 
ies are dimmer than those of the face-on galaxies (which 
shows the combined effect of the magnitude cut excluding 
brighter edge-on galaxies and our size cut which allows in in- 
clined galaxies at observed dimmer magnitudes if they have 
increased observed radii). The intrinsic magnitudes of the 
inclined galaxies must be brighter than their observed mag- 
nitudes, so the distribution of intrinsic magnitudes across 
the inclination range of our sample may be more similar. 
Converting observed axial ratios to inclinations depends 



on assumptions about the intrinsic axial ratio of the galaxy 
(as discussed in Section 2.4.2 above). In order to provide 
empirical corrections, and compare observed data to models, 
we avoid making this conversion (which adds an additional 
source of error), although we provide a prescription for es- 
timating the intrinsic axial ratio from fDeV in Section 2.4.2 
which we will apply here. The cos i distribution of the sam- 
ple is a good sanity check of the level of sample bias. This 
distribution should be flat (i.e. random orientations) if the 
sample is completely unbiased. We show this distribution 
in Figure [5] There is an excess of objects at i ~ 78°, and a 
deficit at i < 25° and i > 84° . However we estimate a typical 
error of ±0.1 in the estimate of cos i (coming from an error of 
±0.1 in the measurement of b/a and ±0.05 in the estimate of 
q) and argue the distribution is reasonably flat within this 
error. For a very conservative interpretation of the results 
shown below the region 25° < i < 70° can be considered 
completely unbiased (roughly 0.04 < log(a/6) < 0.45). 

That this sample ends up being relatively unbiased (as 
shown by Figure 8) can be explained if either the biases 
introduced by the two selection effects (SDSS magnitude 
limit and our size cut) cancel out precisely, or if one of the 
selections dominates and this selection is relatively unbi- 
ased. The first explanation, while possible, seems unlikely. 
In fact, Figures 6 & 7 seem to show that the size selection 
is the dominant selection effect (i.e. most objects not mak- 
ing the SDSS magnitude cut would not have made the size 
cut anyway). This is most easily argued from Figure 6, since 
a shift in r-band magnitude of ~ 0.5-1.0 aligns the size- 
magnitude relations for edge-on and face-on spirals with no 
need for size shifts, and this is in agreement with amount 
of total attenuation estimated in r-band by us and othe r 
authors (see Section 3, and e.g. lUnterboTn fc Rvdenll2008h . 
This then suggests that inclination has little effect on the ob- 
served sizes of spiral galaxies - at least at the radius traced 
by rgo- This is possible, if the galaxi es are mostly opaque 
within this radius in the optical bands. iMasters et al.l (120031 ) 
using 2MASS data, asked if spiral galaxies are transparent in 
NIR bands at 2-3 scale lengths, and conclude this is unlikely 
in the J, and H-bands (but consistent with the data at K- 
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band). The radius, r$o (enclosing 90% of the Petrosian flux) 
is at ~ 4 scale lengths for a pure exponential l|Strauss et al.l 
2002, Figure 1), but it still does not seem inconsistent that 
spirals would still be relatively opaque in the shorter optical 
bands at this radius. It seems therefore that the size cut we 
apply results in a relatively unbiased sample. 



3 EFFECT OF VIEWING ANGLE ON THE 
COLOURS OF SPIRAL GALAXIES 

We observe disk galaxies at a variety of inclinations rang- 
ing from completely face-on (with axial ratios ~ 1, so 
log(a/6) = 0) to completely edge on (with large axial ra- 
tios, e.g. a/b = 10 corresponding to log(a/6) = 1.0). The 
universe is assumed to be isotropic and homogeneous, so we 
expect that the intrinsic properties of galaxies should not 
vary with viewing angle. Any such variation in an unbiased 
sample can therefore be largely interpreted as the effect of an 
increased path length through the disk of the galaxy (with 
complications due to the different distributions of stars and 
dust in a galaxy, which may be averaged out in a large sam- 
ple). The increased path length in the presence of dust will 
result in reddening, and dimming, of the emerging light. 

In this section, we plot the variation of colour with axial 
ratio (as a proxy for inclination) in our sample of GZ spi- 
rals. We use a simple linear parameterization with log(a/i>) 
to quantify the effect of inclination on the observed magni- 
tudes and colours of the galaxies. This is the parameteriza- 
tion that has been historically used in studies of inclination 
dependent effects, presumably arising from the fact that for 
a totally opaque, pure exponential disk the total magnitude 
will change by 2.51og(a/b). We write, 

Xtrue = X obs - 7x log(a/6). (2) 

where X describes the band in question (X = ugriz for 
SDSS bands), and jx for an individual galaxy might rea- 
sonably be expected to depend on things like the galaxy 
mass and /or luminosity, its intrinsic colour, its metallicity, 
its Hubble type, etc. 

In Figure [9] we show the effect of inclination on the 
colour of our spiral galaxies. As expected we find a strong 
trend such that more elongated (or inclined) galaxies are 
reddened. We fit our linear parameterization, as discussed 
above, in the range of log(a/6) < 0.7 and find; 

u-z = 2.267(7) + 0.55(2) log(a/b), 

g-z = 0.969(4) + 0.39(1) log(a/6), 

r-z = 0.418(3) +0.25(1) log(a/6), 

i-z = 0.180(2) +0.12(1) log(a/6). (3) 

The trends of other colours can be found using linear 
combinations of these relationships. Note that the errors 
here are statistical - the systematic error due to the pos- 
sible sample biases discussed in Section 2.5 is likely larger. 

The scale of the y-axis in the different panels of Figure 
[9] has been set by the overall dispersion in the colours of our 
galaxies, and varies from 3 magnitudes in u — z to only 0.7 
magnitudes in i — z. The trend of the reddening due to dust, 
while clearly present in all the observed colours shown here, 
is on average smaller than the overall dispersion seen in this 
sample of visually selected spirals (which is a mixture of 



all Hubble types). It appears that some galaxy colours can 
be massively affected by dust, while others could be affected 
very little at all. This plot illustrates however that the stellar 
populations (mostly through the size of bulge, or old stellar 
component), are just as important in explaining the colors 
of these well-resolved spiral galaxies as the amount of dust 
in their disks. Dust should be considered just one of many 
factors explaining the colours of spirals - however its sys- 
tematic reddening of spiral colours with inclination means 
it must be dealt with if inclined systems are to be included 
in studies of galaxy evolution. 

At log(a/6) > 0.7, the colour trends above are sup- 
pressed by a downturn in the mean colours for the most 
elongated spirals. We interpret this effect as being due to 
the fact that only the latest, and bluest, of the spiral types 
(i.e. Scs and Sds with no bulge component) can have such 
elongated shapes. As discussed below, an early-type spiral 
will reach i — 90° at log(a/6) ~ 0.7 so cannot have any large 
value of a/b. This illustrates the difficulties of using axial ra- 
tio as a proxy for inclination in a sample which includes a 
wide range of spiral galaxy types. The average colours of the 
most inclined galaxies may also be impacted by our se- 
lection which allows in intrinsically smaller galaxies at high 
inclinations, which on average will also be less luminous and 
possibly bluer (through the colour- luminosity relation). It 
could also be observed if the disks become totally opaque at 
these high inclinations and the light is therefore dominated 
by bluer emission from the near side disk. 

Our observed trends of A — z with log(a/6) can provide a 
lower limit to the average value of jx for spiral galaxies (un- 
der the assumption that the internal extinction in z-band is 
negligible). The SDSS z-band should have similar (or slightly 
more) internal extinction than the J-band, which was mea- 
sured as 7 j ~ 0.5 for a sample of vis u ally c lassified spirals 
detected in 2MASS bv lMasters et all ((2003), thus suggest- 
ing 7 Z > 0.5.0 Using this argument, we use the above SDSS 
colour trends to show that, on average (neglecting depen- 
dences on galaxy mass, luminosity, colour, metallicity and 
type), j u > 1-0, 7 9 > 0.9, 7 r > 0.8, ji > 0.6 which we 
suggest should be used to log(a/6) ~ 0.7 after which the 
extinction appears to be on average constant with axial ra- 
tio. This means that a lower limit on the total extinction 
of spiral galaxies from face-on to edge-on, is 0.7, 0.6, 0.5 
and 0.4 magnitudes in ugri bands respectively (based on 
an assumption of at least ~ 0.3 magnitudes of extinction in 
z-band from face-on to edge-on). 

We advise against using uncorrected axial ratios as a 
proxy for inclination in galaxies smaller than rgo = 10" due 
to the effects of seeing which will make the most inclined 
galaxies appear rounder than they should be (as shown in 
Figure 1). This effect will compress the trend of reddening 
with inclination making it appear larger than it really is. 
However we recognize that a correction to the mean colours 
to the face-on value may still be of use in statistical studies 
of galaxy evolution so provide below trends which are for 
"small GZ spirals" (meaning those with r9o < 10" only, or 



5 iMaller et alj (120091 ') argue for much lower extinction in J-band 
(max(7j) = 0.17), but do quote max(7 z ) = 0.53, so assuming 
7 Z > 0.5 is not inconsistent with that work. 
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Figure 9. Shown here are the trends of u — z, g — z, r — z and i — z colours as a function of axial ratio for our sample of well— resolved GZ 
spirals. Greyscale contours are linear in the density of all galaxies. The median in bins of 400 galaxies are overlaid, as are the 25th and 
75th percentiles. The best fit linear relation (to log(a/6) < 0.7) is shown as the straight line; for log(a/fe) > 0.7 we just show a constant 
extension of this fit. Note that the y-axis range, set by the overall dispersion in colour differs dramatically between (u — z) where is is 3 
magnitudes, to (i — z) where is is only 0.7 magnitudes. 



55,653 galaxies in total): 

u-z = 1.980(4) + 1.221(13) log(a/fc) 

g - z = 0.863(3) + 0.692(7) log(o/6) 

r - z = 0.380(2) + 0.388(4) log(o/6) 

i-z = 0.188(1) +0.165(2) log(a/6). 



(4) 



3.1 Dependence on Bulge-Disk Ratio or Spiral 
Type 

In Figure [5] there is evidence for a possible bimodality in 
the trend of the colours with axial ratio. As we will see be- 
low, w hen discussing the photometric model of I Tuffs et al. I 
|2004l ). the bulge-disk ratio (or spiral type) affects the model 
attenuation curves. We also know that the average colours 
of spiral galaxies becomes bluer as they progress along the 
Hubble sequence (from early to late spirals) as their bulges 
become smaller. So an explanation for the apparent split in 
the colour-log(a/&) trend may be t he bulge-disk ratio of t he 
spirals (or perhaps bulge type - see lDrorv fc Fis her 2007). 
As discussed in Section 2.4 above, we can use fnev 



to separate GZ spirals into subsamples based on the size 
of their bulge (or give a crude spiral type classifications). 
Obviously it would be better to use full bulge-disk decom- 
positions to study these effects which would be a natural 
extension of this work. In Figure [lOl we show the trends of 
g — z colour with axial ratio in four bins of foev (see Figures 
IB1IIB3I in Appendix B for other colour trends). In Tabled 
we present the best fit parameters to these relations. While 
for the full sample, we only fit the range Iog(a/6) < 0.7, 
for these subsamples (split by /dev) we fit the full range 
of log(a/&) as the flattening in the median colours seen in 
Figure [9] is less pronounced here - perhaps because of the 
split into rough spiral types. 

Figure [lOl shows that as we move from late-type spirals 
(with small bulges) to early-type spirals (with large bulges) 
the average face-on colours (Co in Table [2} reddens sig- 
nificantly. In fact the average face-on colours of the "very 
bulgy" subset are redder than the average edge-on colours of 
the "pure disk" subset (by 0.3, 0.1, 0.01 and 0.002 mags in 
(u — z), (g — z), (r — z) and (i — z) respectively) - although the 
spread of the colours of even face-on "disky" spirals is large 
enough to almost encompass the spread of face-on "bulgy" 
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Figure 10. The trend of g — z colour for GZ spirals separated using the foe.v parameter as indicated in the panels. As in Figure [5] 
the greyscale contours represent the locus of galaxies, while overlaid are the median and interquartile range of data binned in log(a/6) 
(bin sizes are (1) f DeV < 0.1; 163, (2) f DeV < 0.5; 352, (3) f Dt , v > 0.5; 134, (4) f Df , v > 0.9; 50 - designed to make ~ 50 bins for each 
subset) 



Table 2. Best Fit linear Trends to (A — z) = Co + 7log(o/6). Aq is the total reddening from face-on to edge-on. 



Sample 






(u - z) 






to-*) 






(r - z) 






(i-z) 








Co 


7 


A 


Co 


7 


A 


Co 


7 


A 


Co 


7 


A 


All big spirals 


24276 


2.267(7) 


0.55(2) 


0.39 


0.969(4) 


0.39(1) 


0.27 


0.418(3) 


0.25(1) 


0.18 


0.180(2) 


0.12(1) 


0.08 


"Very bulgy" 


2519 


2.81(1) 


0.62(4) 


0.43 


1.257(7) 


0.39(2) 


0.27 


0.575(4) 


0.22(1) 


0.15 


0.253(3) 


0.11(1) 


0.07 


"Bulgy" 


6680 


2.68(1) 


0.74(3) 


0.52 


1.198(5) 


0.48(1) 


0.34 


0.547(3) 


0.28(1) 


0.20 


0.238(3) 


0.15(1) 


0.11 


"Disky" 


17596 


2.030(8) 


0.67(2) 


0.67 


0.842(4) 


0.45(1) 


0.45 


0.349(3) 


0.28(1) 


0.28 


0.146(2) 


0.14(1) 


0.14 


"Pure disk" 


8168 


1.91(1) 


0.57(2) 


0.57 


0.766(7) 


0.40(1) 


0.40 


0.303(4) 


0.26(1) 


0.26 


0.121(4) 


0.13(1) 


0.12 



spiral colours. As discussed above, dust reddening appears 
again to be about equal in predicting the colors of galaxies 
as the bulge size and the stellar mix of the galaxy (however 
the difference is its systematic trend with viewing angle). 

We observe that the slope of the reddening with axial 
ratio is steepest for the "bulgy" subsample. If we assume 
that a "disky" spiral (i.e. one with fDeV < 0.5) has an 
intrinsic axial ratio of q ~ 0.1, while a "bulgy" spiral has 
q ~ 0.2 (see Table 1 for a justification of these choices), then 
the total reddening from face-on to edge-on (noted as Ao in 
Table [2} is given by 7 for the "disky" (late-type) spirals 



and by O.77 for the "bulgy" (early-type) spirals. Using this 
assumption, we argue that the total reddening from face-on 
to edge-on is always larger in the "disky" spirals than in 
the "bulgy" ones which then indicates a larger role for dust 
reddening on the average colours of the "disky" spirals than 
in the "bulgy" spirals. 



We test this data against the photometric models of 
iTuffs et al.l (|2004h in Section 4 below. 
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3.2 Dependence on Galaxy Luminosity 

Luminosity dependence of the amount of du st attenuation in 
spiral gala xies was studied in NI R bands in iGiovanelli et al.l 
il 19951) and lMasters et~liIl (l2003h and in the optical/NIR in 
iTullv et all (|l998l ). iMaller et al.1 (120091 ) have also studied 
the impact of luminosity and face-on colour on reddening 
trends of late-types in SDSS bands. In the previous section 
we showed differences in the trends of reddening with incli- 
nation for GZ spirals with different bulges sizes (as measured 
with JdeV from SDSS). It has been shown many times that 
different types of spirals have different luminosity functions 
{e.g. work using 1500 SPSS galaxies with visual classifica- 
tion, iNakamura et al.l 120031 ). so some, or all of this trend 
could be due to the different luminosity ranges of spirals 
with different bulge sizes. 

In Figure [TT1 we show the observed (i.e. not corrected 
for internal extinction) colour magnitude (CM) diagram of 
all GZ spirals (p sp ira,i > 0.8, 0.01 < z < 0.09) in four quartile 
bins of axial ratio. Each plot has ~ 20000 galaxies on it. The 
galaxies are colour coded by fDeV using: red: /deV > 0.9; 
yellow: 0.5 < /o e v > 0.9; green: 0.1 < /o e v > 0.5; and 
blue: fDeV < 0.1. The solid line indicates the blue edge of 
the red sequence of GZ early-types at z = 0. The "face-on" 
GZ spirals are mostly to the blue side of this line, but in the 
most edge-on quartile many galaxies are reddened across 
it. We note that most spirals whose intrinsic colours place 
them in the so-called "green valley" are bulge-dominated, 
but there is significant contamination to this part of the 
CM diagram by dust reddened inclined spirals. 

The dashed line in each panel of Figure QT] is a fit to 
the blue cloud in the face-on GZ spirals (log(a/6) < 0.14) 
- with the vertical dashed line indicating the 99.5th per- 
centile magnitude (i.e. only 5% of spirals have r-band mag- 
nitudes brighter than that line). The dotted lines in the 3 
panels showing more inclined spirals are the same quantities 
for those inclined galaxies and illustrate the change in the 
shape of the CM diagram of spirals as they become more 
inclined. In the most inclined quartile the 95th percentile 
r-band magnitude is 0.5 mag fainter than in the face-on 
spirals, quite consistent with the estimate of the trend of 
r-band magnitude with axial ratio we derived in Section 3.1 
above. 

We derive the face-on r-band absolute magnitudes of 
the well resolved GZ spirals using a correction of AM r = 
0.81og(a/6) to log(a/6) = 0.7, and AM = 0.56 after (as sug- 
gested in Section 3) then fit linear relations to the colour ver- 
sus axial ratio in subset of the GZ spirals separated by both 
bulge size (as measured by fDev) and absolute magnitude. 
We plot the slopes of these relations for the (r— z) colour ver- 
sus r-band magnitude in Figure [T2l Interestingly the slope of 
the trend of reddening with axial ratio shows hints of a peak 
in all spiral type subsets at around M r ~ —21.5 where it is 
roughly twice as large as at the lowest luminosities. At both 
dimmer and brighter magnitudes the reddening trend seems 
to decrease (note that in the subset of the most bulge dom- 
inated spirals there are very few galaxies with M r > —21). 
We comment that while we have selected the galaxies in 
these magnitude bins after an estimated dust correction has 
been applied, the implicit magnitude limit of the original se- 
lection may be causing some bias here (i.e. causing edge-on 
galaxies which would make it into the sample if seen face-on 



to be dimmed below the selection limit) and would act in 
a way to cause an underestimate of the slope of reddening 
(since brighter galaxies will on average be redder). However 
the blue cloud shows only a mild trend of colour with lumi- 
nosity, so we argue that this effect is unlikely to account for 
all of the trend with luminosity we observe here. However 
the bias may be larger at the faint end of the sample, so the 
downturn in reddening at in low luminosity spirals should 
be interpreted carefully. 

At a given absolute magnitude the slope of the redden- 
ing is significantly shallower (0.1-0.2 mag/dex) in the "very 
bulgy" subset /dcV > 0.9 than in more disk dominated spi- 
rals. This means that a single correction to magnitudes and 
colours for dust will be an overestimate of the effect in both 
very dim and very bright (especially bulge dominated) spi- 
rals and an underestimate especially for disk dominated spi- 
rals with M r ~ —21.5. Interpreting these slopes in terms of 
dust content is not trivial as illustrated in Section 4. The 
slope of the predicted tr end of (r — z) colour from the mod- 
els of lTuffs et al~l (2004) is degenerate with bulge-disk ratio 
and face-on central opacity. However in pure disk galaxies 
which our /zjev < 0.1 sample should approximate (shown 
in blue in Figure [12} an increase in the slope of the trend of 
reddening can be interpreted as an increase in dust content 
in the galaxies. 

An increase in dust content is expected for more lumi- 
nous galaxies which have had more star formation over cos- 
mic time and which also are physically larger and so their 
path lengths increase. A decrease of the trend for the most 
luminous galaxies has not been observed before, perhaps 
because previous studies have focussed on later type spirals 
which are dominated by the lower luminosities. However it 
has been observed that lower luminosity galaxies have more 
active recent star for mation than high mass galaxies even 
amon gst spirals (eg. Ybund 19991 ). Dust is destroyed over 
time |Draine fc Salpetei 19791 ). so these hints in the down- 
turn in the slope of reddening in the most luminous galaxies 
may be related to their lower levels of recent star formation. 

Clearly the total reddening of a spiral galaxy with in- 
clination is a complex mixture of it's dust content (via its 
star formation history), physical size and dust geometry, 
all of which may depend on other global quantities such 
as luminosity, or morphological type. In order to under- 
stand dust attenuation we therefore argue a two pronged 
approach is needed. Individual modeling of galaxies using 
multi-wavelength data (from UV to far-IR) is needed to shed 
light on the details of the dependence of dust on galaxy 
properties; while the kind of statistical approach used in 
this paper is needed to look at global trends in large sam- 
ples of galaxies, and to provide empirical corrections of use 
to remove, at least to first order, the bias introduced into 
galaxy surveys from the inclination dependent attenuation 
and reddening of spiral galaxies. 

3.3 Galaxy Zoo Red Spirals 

In lBamford et al.1 12009), we discuss the existence of a large 
population of "red spirals" detected in the GZ sample, 
i.e., galaxies with colours consistent with elliptical galax- 
ies ((« — r) colour gr e ater t han the dividing line in Fig- 
ure 7c of lBaldrv etaU l|2006h ) but with clear evidence of a 
disk and/or spiral arms. We revisit these "red spirals" here 
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Figure 11. Colour- magnitude (M r vs. (g — r)) for Galaxy Zoo spirals in four quartile bins of axial ratio (as indicated). Each plot has 
~ 20000 galaxies on it. Neither colour nor magnitude is corrected for internal extinction. The galaxies are colour coded by fDeV using: 
red: foeV > 0-9; yellow: 0.5 < fDeV > 0-9; green: 0.1 < foeV > 0-5; an <3 blue: foeV < 0.1. The solid line indicates the blue edge of the 
red sequence of GZ ellipticals. The dashed line in each panel is a fit to the blue cloud in the most face-on quarter of Galaxy Zoo spirals 
(log(o/6) < 0.14) - with the vertical dashed line indicating the 99.5th percentile magnitude (i.e. only 5% of GZ spirals in this axial ratio 
bin have r-band magnitudes brighter than that. The dotted lines in the 3 panels showing more inclined spirals are the same quantities 
for those galaxies and illustrate the change in the shape of the CM diagram of spirals as they become more inclined. 



mainly to note that the dividing line between "red" and 
"blue" changes as a function of the derived stellar mass of 
the galaxy, the measurement of which can depend on spi- 
ral galaxy inclination. The presence of dust both increases 
the mass-light ratio (redder colours imply higher M/L) and 
decreases the luminosity (due to extinct ion) so a trade off 
exists as discussed by other authors (e.g. iDriver et ai1l2008l ; 
iMaller et aHl2009h . In lSkibba et al.1 i|2009h . we presented the 
clustering of the red spirals (now selected from the M r ver- 
sus (g — r) colour-magnitude and examined the expected 
contamination of this sample by reddened edge-on spirals. 

We emphasize that we still find a population of inter- 
esting "red spirals". For example, we see face-on GZ spi- 
rals that are still redder than the normal colour dividing 
line between spirals and ellipticals (eg. Figure [TTJ, and even 
the reddest edge-on spirals must be unusually passive un- 
less they have signifi cant amounts of dust. Furthermore, 
iBamford et al.l l|2009l ) demonstrated that the face-on red 
spirals have similar environmental dependencies as all red 



spirals (including the dust reddened spirals) . [Masters et al. I 
(20101) study in more details these intrinsically red spirals. 



4 PHOTOMETRIC MODELS FOR DUST 
ATTENUATION IN SPIRAL GALAXIES 

We next compare models for dust attenuation with our ob- 
served relations. There has been significant progress in mod- 
eling the attenuation of the total light from spiral galaxies 
as a function of i nclination, and mode ls continue to grow in 
complexity. Both lFerrara et al.1 (|l999l ) and Tuffs et al. 2004 
(hereafter T04) provide models which allow the attenuation 
to be calculated for a spiral galaxy at any viewing angle and 
with a given bulge-to-disk ratio. For a history of dust mod- 
eling in spira l galax ies see the introduction of either T04 or 
IPierini et al.l l|2004l ). A more recent discussio n of the T04 
model can be found in iPopescu fc Tuffsj i|2009l ). 

Here we compare our data to the T04 model. In this 
model, a galaxy can be constructed with any bulge-disk ra- 
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Figure 12. Slope of the trend of (r — z) reddening with axial ra- 
tio plotted against the r-band absolute magnitude (with a basic 
correction for dust as described in the text). The different lines 
show subsets divided by foeV as used elsewhere in this work go- 
ing from blue-green— orange— red (solid-dashed-dot-dash-dot-dot- 
dot-dashed; symbols square-diamond-triangle-circle) from most 
disk dominated to most bulge dominated. 



tio. The attenuation of each of three components (a bulge, 
disk and thin disk) is provided as a functional fit with in- 
clination at central optical depths in B-band of tb = 0.1, 
0.3, 0.5, TO, 2.0, 4.0 or 8.0. The different components are 
then combined in the appropriate fractions to give the total 
attenuation. The output of the T04 model is published in 
various UV bands, as well as BVIJK. We linearly interpolate 
these results to the central wavelengths of the SDSS ugriz 
bands (as recommended by Tuffs, 2009, priv. comm.). 

When attenuation curves are considered {i.e. attenu- 
ation as a function of inclination) the relative geometry of 
the dust and stars are just as important as the wavelength 
dependence of the extinction (which is "Milky Way like" in 
T04). Details of the exact distributions and relative sizes 
of the diffuse stellar light and dust can be found in T04. 
For our application it is important to know that the relative 
geometry is fixed based on optical/NIR imaging of a small 
sample of edge-on galaxies (jXilouris et al.l ll999). and that 
also the relative face-on B-band optical depth of the thick 
and thin disks is fixed at a r atio of 0.387 as mea sured for 
the edge-on galaxy NGC 891 (|Popescu et alJl200Ch . 

The free parameters are then the total face-on B-band 
opacity (r = r d i s k + Tthindisk), how "clumpy" the clumpy 
dust component is (expressed by the dumpiness factor, F, 
for which a suggested value of F = 0.22 is given in T04 
based on observations of NGC 891), and the inclination of 
the galaxy. The observed attenuation also depends on the 
observed bulge-disk ratio of the galaxy. 

One final ingredient which will be important in compar- 
ing the T04 model to our data is the dependence of intrinsic 
axial ratio on the intrinsic bulge-disk ratio. Only in the case 
of an infinitely thin disk is the inclination related to the ax- 
ial ratio by cos(i) = b/a. In normal disks with an intrinsic 
axial ratio, q, the inclination is given by Equation [T] and 
we expect that q will increase monotonically with increasing 
bulge-disk ratio. 
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Figure 13. The average reddening A — z from face-on (i = Odeg) 
to edge-on (i = 90 deg) from the models of T04. Top panel: the 
results for a pure disk model (B/D = 0, solid line) and a large 
bulge-disk ratio (B/D = 6; dot-dot-dot-dashcd line) for central 
B-band face-on opacities of 0.1, 0.3, 0.5, 1, 2, 4 and 8 (from lower 
to upper). Middle Panel: the results for varying the bulge-disk 
ratio from pure disk (B/D = 0) to large bulge-disk ratios (up to 
B/D = 6) when rg = 2. Lower Panel: the result of changing 
the dumpiness factor, F with tb = 2 and B/T = 0.3 (Sb like 
spirals). Overlaid with symbols are the results of our fits to the 
"very bulgy" , "bulgy" , "disky" and "pure disk" subsets of the GZ 
spirals (circles, squares, triangles and stars respectively; see Table 
0. 



4.1 Direct Comparison 

In Figure [13] we plot lines showing the dependence of the 
total reddening from edge-on to face-on with wavelength for 
the free parameters in the T04 model. The top panel is an 
illustration of the impact of changing the central B-band 
face-on opacity in the case of a pure disk model (solid line) 
and a galaxy with a large bulge-disk ratio (B/D — 6, dot- 
dash line). In the middle panel the impact of changing the 
bulge-disk ratio from a pure disk galaxy through to large 
bulge-disk ratios for a model with t% = 2 is shown. In the 
bottom panel we show the impact on (u — z) reddening of 
the full range of possible values for the dumpiness factor F. 
This factor only affects the u-band in the SDSS filter. 

Also plotted in Figure [T3] are fits to the observed to- 
tal attenuations from Section [3] While the model can easily 
predict the attenuation in the optical colours, it is clear that 
the edge-on attenuation curve significantly over-predicts the 
observed edge-on u-band attenuation for reasonable values 
of the central B-band opacity. The u-band is the only SDSS 
band in which the thin disk contributes to the attenuation, 
which causes the attenuation to increase significantly only 
at inclinations where the thin disk is close to being viewed 
edge on. The exact wavelength transition between where the 
thin and thick disks are important is not well known. The 
T04 model can be modified to make the it-band be domi- 
nated by the thick disk (by multiplying the B-band output 
of the model by 1.267) and this provides a better fit to the 
GZ spirals (Tuffs & Popescu, 2009, priv. comm.). 

We show Figure[l4]as an example of the effect of varying 
model parameters on the curves of optical reddening with 
axial ratio. We show this for a range of central B-band opaci- 
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Figure 14. The (r — z) reddening as a function of observed axial 
ratio from the models of T04. Top panel: the results for a pure 
disk model (B/D = 0) for central B-band face-on opacities of 
0.1, 0.3, 0.5, 1, 2, 4 and 8 (from lower to upper). Lower panel: the 
results for t b = 2 and B/T values of 0, 0.25, 0.5, 0.75 and 1.0. 

ties (top panel) and for a range of bulge-disk ratios at tb = 2 
(lower panel). It is interesting to note that while increasing 
the central opacity (or the total amount of dust) increases 
the total amount of attenuation significantly, the slope of the 
curve of attenuation with axial ratio does not change much. 
Unfortunately this is the main observable we are sensitive to 
(not knowing a priori the expected intrinsic colours of spi- 
ral galaxies), which (if the models are correct) will limit the 
information on tb which can be obtained from studying the 
inclination dependent optical attenuation of spiral galaxies 
(at least via colours). The effect of an increasing bulge-disk 
ratio increases the slope of the curve since at larger B/D 
more inclined galaxies have smaller axial ratios. Also strik- 
ing is the downturn in attenuation at large inclinations when 
bulges are present (as discussed in T04). 

The amount of face-on reddening from the T04 models 
is only a fraction of a magnitude even for tb = 8, showing 
that (even given the large trends of observed colour with 
axial ratio), dust reddening of face-on galaxies cannot be 
significant. The range in face-on colours predicted by T04 
from central opacities varying from tb = 0-8 is also much 
smaller than the observed range of face-on colours of galax- 
ies {e.g. Figure [3|. This illustrates that (if this model is 
correct) the colour of a face-on spiral galaxy is set by the 
intrinsic colour of its stellar population and not by the red- 
dening effects of dust. This therefore shows (again if the 
T04 model is c orrect) that the face -on red spirals observed 
in Galaxy Zoo (|Masters et al. ifeoiCf ) must have intrinsically 
red colours from old stellar populations. 

4.2 Monte Carlo Realisations 

To produce a fair comparison of our data with the models, 
we want to fully sample the intrinsic scatter in the proper- 
ties of our GZ spirals. Therefore, we produce Monte Carlo 
realisations of the model attenuation curves in the SDSS 
bands based on reasonable assumptions about the galaxies. 
We use the correlations found in Section \2. 41 between JdsV, 
bulge-disk ratio, intrinsic axial ratio and face-on colour to 



construct Monte Carlo estimates of the distribution of ob- 
served colours from the T04 models. 
The Monte Carlo ingredients are: 

1. tb - assumed to be Gaussian with a mean tailored to 
fit the data and standard deviation of 1. 

2. F - has no affect on attenuation curves, so set to 
F=0.22 

3. B/T - modelled using the observed distribution of 
fDeV and the relation between foev and B/T we find in 
Section T2.4I (adding in a scatter of 0.1 in B/T). 

4. log(a/6) - from Eqn[T] assuming random inclinations 
and intrinsic axial ratios which depend on fDeV as found in 
Section |2~41 2 

5. Face-on colours - input as a function of fDeV based on 
observed colours of face-on GZ spirals (Figure 3). 

We try various input values for tb- For low values of 
tb (less than tb ~ 1.5 ± 1) the model is able to predict the 
observed trend of (u — z) colours, but significantly under- 
predicts the optical colours of the most inclined galaxies. For 
any value of tb larger than about 1.5, the (u — z) colour of 
the most inclined galaxies is over-predicted while the optical 
colours of the most inclined galaxies are under-predicted. 
(The over-prediction in it-band may be improved using the 
thick disk modification to the B-band attenuation in the T04 
model which is discussed above.) 

We do not reproduce the shape of the scatter in Figure 
2, which shows hints of a bimodal pattern, and has a long 
trail to the blue side. The scatter in the Monte Carlo is dom- 
inated by our input assumption of a Gaussian distribution 
of face-on colours, so this shows a more complicated distri- 
bution of face-on colours is needed. The interquartile ranges 
matches well at low inclinations (which is expected as this is 
an input based on the data) but at large inclinations - where 
the effect of the model start to dominate - the interquartile 
range it is found to be too large in the M-band colour and too 
small in the optical colours. This suggests that our Monte 
Carlo using the T04 model predicts a larger range of UV 
reddening and a smaller range of optical reddening than is 
observed. 

Figure [15] shows our model Monte Carlo realisations 
for the reddening as a function of inclination with an input 
central face-on B-band opacity of tb = 4 ± 1. We pick this 
value of tb i n par t based on the tb = 3.8 ± 0.7 fou nd by 
iDriver et all (|2007h using the same lTuffs et al. I (|2004T ) mod- 
els and looking at the inclination dependence of the B-band 
luminosity function, but it does also appear to describe the 
observed trends well. We note here that comparisons of our 
data with this model are therefore implicitly also compar- 
is ons between the the trends we present and those presented 
in IDriver et all (|2008f l which are based on the T04 model 
with this value of tb- We overplot (in orange) the binned 
real data and linear fit to that data. In purple, we show the 
same quantities but for the Monte Carlo realisations. 

We split the sample by fDeV as done in Section 3.1 (see 
Figures [T0l and [BlHB3l for the data). Here we just explore the 
two extreme subsets, the "pure disks" (fDeV < 0.1) and the 
"very bulgy" (/dcV > 0.9) spirals (see Figures [Cll and [C2)l . 
This split removes worries about the conversion between in- 
clination and observed axial ratio, as all galaxies in these 
two extremes of bulge size should have similar observed ax- 
ial ratios at a given inclination. 
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In summary we find that the model reproduces the ob- 
served trends surprisingly well, considering the number of 
fixed geometrical parameters (eg. the relative size of the 
stellar and dust disks, bulge shape, etc.) and agrees with 
the value of tb = 3 .8 ± 0.7 (central optical depth) found by 
iDriver et al. (2007) for the same models. However, the mod- 
els overpredict the amount of u-band attenuation in edge-on 
galaxies while underpredicting in gri-bands for any reason- 
able choice of central opacity and bulge-disk ratio (but see 
previous comments about modifying the published model to 
better match our it-band data). This could point to an inad- 
equacy of the Milky Way extinction law when applied to ex- 
ternal galaxies, but more likely indicates that a wider range 
of model geometries need to be allowed for. The relative ge- 
ometry of the dust and stars is as important as total dust 
content when considering the integrated colours of galaxies, 
an d may differ quit e subst antially among spirals. For exam- 
ple iHolwerda etafl (120091 1 study the dust distribution in a 
backlit galaxy observed with ACS, and find it has a very ex- 
tended dust distribution, clearly illustrat ing the limi t ations 
of using a fixed geometry for all galaxies. IWitt et al.l (|l992T l 
show how very dusty galaxies can be quite blue if the stars 
extend well past the dust so that there is a blue contribu- 
tion from the outer disk. The same effect might be observed 
if the dust is very patchy. The observed light will become 
progressively dominated by emission from stars in the clear 
regions as the wavelength shortens. This could easily explain 
the failure of the T04 model in the SDSS it-band. The in- 
creased scatter in the reddening might also be related to the 
patchy distribution of dust in spiral arms and the expected 
variation of the relative orientation of spiral arms with the 
major axis of the inclined spiral. 



5 SUMMARY AND CONCLUSIONS 

We present here a careful study of 24,276 "well-resolved" 
(significantly larger than the average SDSS seeing) spiral 
galaxies selected from the visual classifications of the Galaxy 
Zoo project. This sample avoids the effect of seeing on the 
axial ratios of these galaxies, thus allowing us to use them 
as a reliable proxy for inclination (see Figure 1). We then 
study the properties of these spirals as a function of inclina- 
tion to infer the effects of galactic dust on these properties. 
We also compare the data with the latest models for dust 
attentuation. Our major results and conclusions are: 

• We find the SDSS parameter fDeV (fraction of the light 
profile best fit by a de Vaucouleurs model) is a good can- 
didate for estimating the bulge size (and thus spiral type) 
for our galaxies, and appears to be better than other mor- 
phological proxies commonly used in the literature (e.g. con- 
centration has a significant dependence on axial ratio with 
more inclined spiral galaxies being more concentrated). By 
matching to a sample of MGC galaxies (with known B/T 
ratios) we find that JdhV correlates with the bulge size with 
a linear relationship of fDeV = 0.1f(3) + 1.7(2)£>/T. As a 
confirmation of this result, we see that the face-on colours 
of GZ spirals show a clear trend with /nev such that galax- 
ies with larger fDeV are redder as expected for spirals with 
larger bulges (see Figure [2}. 

• We observe a clear correlation between the colours of 
our spirals and the inclination of the disk, i.e., those with 



larger axial ratios (more inclined) are on average redder than 
those which are viewed face-on (with small axial ratios). We 
have fit these relations and provide them in Section 3.1. We 
then use these fits to estimate the total average extinction 
of spiral galaxies from face-on to edge-on, which is 0.7, 0.6, 
0.5, 0.4 magnitudes in ugri bands respectively (assuming 0.3 
magnitudes of extinction in z-band). 

• Using fDeV, we have split the spiral sample into "bulgy" 
and "disky" spirals (or early- and late-type spirals). We see 
that the average face-on colours of the "bulgy" spirals is red- 
der than the average edge-on colours of the "disky" subsam- 
ple (by as much at 0.3 mags in (u — z)). Likewise, the scatter 
in colors seen in Figures 3 and 6 is slightly greater than the 
trends with inclination or increased reddening. Therefore we 
show that the colours of a spiral galaxy are set about equally 
by the bulge-to-disk ratio, or spiral type (via the stellar pop- 
ulations in the disk and bulge) and reddening as a function 
of inclination, due to dust. However we comment that dust 
effects are systematic with inclination and therefore need to 
be accounted for to make unbiased samples. 

• We explore the luminosity dependence of the redden- 
ing and show that the slope of the trend peaks in all types 
of spirals (as split by fDev) at around M r ~ —21.5. An 
incre ase in the slope with luminosity has been observed be- 
fore dGiovanelli et al.lll995l; iTullv et al. 1998:;; iMasters et al.l 
2003). An increase in dust content is expected for more lu- 
minous galaxies which have had more star formation over 
cosmic time and which also are physically larger and so their 
path lengths increase. A decrease of the trend for the most 
luminous has not been observed before, and is mostly likely 
caused by the lower levels of recent star formation which are 
observed in the most massive spirals (since dust is destroyed 
over time). 

• We compare our data with state-of-the-art dust atten- 
uation models from T04. We construct a Monte Carlo real- 
isation of the GZ spiral samples using reasonable estimates 
for the input range of inclinations (assumed to be random), 
bulge-disk ratios, intrinsic axial ratios and face-on colours 
(taken from the observed distribution of fDev)- We find that 
the model reproduces the observed trends surprisingly well, 
considering the number of fixed geometrical parameters (eg. 
the relative size of the stellar and dust disks, bulge shape, 
etc.) and agrees with the value of tb = 3 . 8 ±0. 7 (central op- 
tical depth) found bv lDriver et al.1 (|2007l . 120081 ) for the same 
models. However, the models overpredict the amount of u- 
band attenuation in edge-on galaxies while underpredicting 
in gri-bands for any reasonable choice of central opacity 
and bulge-disk ratio. They also are not able to predict an 
increase in the range of the observed colours of galaxies with 
inclination. This could point to an inadequacy of the Milky 
Way extinction law when applied to external galaxies, but 
more likely indicates that a wider range of model geometries 
need to be allowed for. In fact modifying the T04 model to 
make the u-band emission come from the thick disk instead 
of the thin disk does appear to fit our data better (Tuffs & 
Popescu, 2009, priv. comm.) 

We end by emphazising the effect of inclination (and 
dust) on present and future large galaxy surveys. As shown, 
the presence of dust in a galaxy both reddens and dims the 
emitted light, and this effect clearly increases for inclined 
spirals. This impacts any measurements based on the mag- 
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Figure 15. Monte Carlo realisation of the reddening in SDSS colours of 24376 spiral galaxies based on the T04 model with Tg = 4 ± 1. 
This figure should be compared to Figure [9] which shows the same thing for the real data. We over-plot in orange the median and 
interquartile range of the binned real data as well as a linear fit to the real data. In purple is shown the same thing for the Monte Carlo 
generated galaxies. 



nitudes and colors of galaxies, like stellar mass estimates, 
but more crucially can cause systematic incompletenesses in 
galaxy survey. In the era of high precision cluste ring mea- 
surements from galaxies (e.g. iPercival et ai1l2009l ). such ef- 
fects could become important to ensure a fair sample of 
galaxies is obtained, e.g., both volume and mass-limited 
samples will miss faint (low mass) inclined spirals which 
fall below the selection criteria because of dust extinction. 
This could induce small bias in the measurements of scale- 
dependent biasing of galaxies as the mass of spirals is corre- 
lated with environment (and preferentially in the field com- 
pared to ellipticals). 
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APPENDIX A: CONCENTRATION AND 
LIGHT PROFILE SHAPE VERSUS GALAXY 
ZOO MORPHOLOGIES 

Without access to Galaxy Zoo, or other visual classifications 
for the huge number of galaxies in SDSS it has become com- 
mon practice to use concentration or other structural pa- 
rameters (for example fDeV which describes the fraction of 
the light profile best fit by a de Vaucouleurs profile as op- 
posed to an exponential profile) as a way to separate early 
and late type galaxies. In this Appendix we show the re- 
lation between these two structural parameters and visual 
classifications from Galaxy Zoo as a guide to ease compar- 
isons between samples split using concentration or /o e v and 
Galaxy Zoo morphologies. 

The first large comparison of fDeV and concentra- 
tion versus morpholo gy in SDSS galaxies can be found in 
IStrateva et alj l|200lf ). There a sample of 287 visual classi- 
fied gal axies and 500 spec t rally classified galaxies were con- 
sidered. IShimasaku et all (|200ll ) performed a similar study 
using a sample of 456 bright galaxies. Both of these papers 
prefer the use of concentr ation (c = rgp/r ^ o) as a morpholog- 
ical separator over juev- Strateva et al. (2001) recommend 
c ~ 2.6 as a separator, while Shimasa ku et alj (|200ll ) pre- 
fer c ~ 3. Both papers warn over the use of fDeV which in 
the brightest galaxies is dominated by the inner light pro- 
file shape. This means that large nearby spirals with bulges 
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Figure Al. The percentage of the "clean" GZ sample (i.e. that 
with either p 3p i ra i or p e u > 0.8) falling into either the spiral 
(blue solid line) or elliptical (red dashed line) classification as a 
function of Top: foeV (the fraction of the light profile fit with 
a de Vaucouleurs profile as opposed to an exponential disk) and 
Bottom: the concentration of the light. 



can have much larger values of fpev than if the fit were 
averaged over the whole galaxy. Never-the-less fpev has be- 
come a common parameter used as a proxy for mo rphology 
{e.g. lUnterborn fe Rvde dliooj ; iMaller et alj|2009h perhaps 
because of the recognition that inclin ed spirals are more con- 
centrated than those viewed face-on (M ailer et alj |2009. and 
see Section T2.4I above). 

Figure [All shows the fraction of the "clean" GZ sample 
(i.e.. those galaxies with very reliable spiral or elliptical clas- 
sifcations; p sp irai > 0.8 or p e \\ > 0.8) which are spirals (blue 
line) or ellipticals (red lines) as a function of fp e v or con- 
centration. We see that while it is true that most ellipticals 
have fp e y ~ 1, many GZ spirals also have quite large val- 
ues of fDeV (as is also evident in Figures [2] and U in Section 
2.4 above). The lower panel shows the same fractions ver- 
sus concentration and shows that the sample is dominated 
by spirals for c < 2.7 — 2.8 and ellipticals at larger con- 
centrations, but that there are still significant numbers of 
galaxies with very reliable spiral classifications from Galaxy 
Zoo (recall these all have p ap irai > 0.8) with relatively large 
concentrations. 

Figure [A"2l shows the cumulative fractions of samples of 
early types or spirals selected by limits on /o e v or concen- 
tration. Again red and blue lines show the fractions for the 
early types and spirals respectively. Fo r example for an fp e y 
separator of 0.5 (as recommended by IStrateva et alj 1200 ll ) 
this Figure shows that roughly 90% of the early types are 
found, and 65% of the spirals. The dotted lines show the 
percent early-type or spiral contamination in the samples. 
For this same example 45% of the "early types" found by 
fDeV > 0.5 are GZ spirals, while 5% of the "spirals" found 
by fpev < 0-5 are GZ early types. In the terminology of 
IStrateva et al.l (|200ll ) this shows the completeness and 100% 
minus the reliability of the samples, and at given values of 
fDeV or concentration can be compared directly with the 
values in their Tables 2 and 3. 

We find that fpev is never a very reliable indicator of 
early types. The highest reliability which can be found is 



Figure A2. Cumulative fractions for the fraction of GZ spirals 
(blue solid line) or ellipticals (red dashed line) found using lim- 
its on Top: foeV or Bottom: concentration. For example, the 
upper plot shows that fp e y < 0.5 will find about 70% of all 
the spirals, while fpeV > 0.5 selects about 90% of the ellipti- 
cals. The orange dotted lines show the elliptical contamination to 
the "spiral" samples, while the purple dot-dashed lines shows the 
spiral contamination to the "elliptical" sample. For example the 
upper plot shows that for fp e V < 0.5, about 10% of the selected 
galaxies are GZ ellipticals, and for fp e y > 0.5 about 40% of the 
selected sample are GZ spirals. 



for fpeV > 1 in which roughly 70% of the GZ early types 
are found, and the contamination by spiral types is 25%. 
Furthermore GZ spiral galaxies with large val ues of fp e y are 
not al ways bright nearby galaxies as found bv lStrateva et all 
l|200ll ). We find GZ spirals with large fpev at the faintest 
magnitudes and largest redshifts in the GZ "clean" sample 
(r ~ 17.7 or g ~ 19, and z ~ 0.25). On the other hand 
fpev can be used fairly safely to identify spiral galaxies. 
Any cut on fpev < 0.6 or so results in a sample which is 
95% GZ spirals and depending on the exact cut will have up 
to 70% of the spirals. The most stringent cuts will miss large 
fractions of the spiral population, and (based on Figure HJ 
will preferentially select late r and later type spirals (i.e. with 
smaller bulges). For example llJnterborn fc R vdcn (200||) use 
fpev < 0.1 to select only the very latest type spirals, but 
we argue that even fpev < 0.5 will miss the earliest spirals 
with the largest bulges. 

Selecting early types by concentr ation is also shown to 
be tricky. The c > 2.6 suggested by IStrateva et alj (|200ll ) 
finds roughly 90% of the GZ early types, but also results in 
a high contamination by GZ spirals (roughly 45%). Many of 
these very concentrated spirals are likely to be edge-on and 
reddened, so an additional colour cut is unlikely to remove 
all of them without also a limit on the axial ratio (as also 
I2009T ). 



su ggested by Maile r et al 



bv lShimasaku et al.l (|200ll ) finds 55% of the early types, but 
results in a sample which is 25% spiral. We find again that 
using concentration cuts to find spirals is more reliable, c < 3 
finds 95% of the GZ spirals with 20% contamination by early 
types, while c < 2.6 finds 75% of the spirals with a 10% 
contamination. 

With the large samples of galaxies now available (and 
planned in the future) there will always be a place for divid- 
ing galaxies in late/early subsets using structural parame- 
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ters. In this Appendix we hope to have given greater insight 
into the levels of completeness and reliability which can be 
expected when using JdbV and concentration (r$o/r$o) from 
SDSS to do this. With Galaxy Zoo 2 data we will be able 
to further study this - looking at the types of spirals which 
end up contaminating the early-type subsets. 



APPENDIX B: TRENDS OF SDSS COLOURS: 
EXTRA FIGURES 

We show here the trends of the u — z, r — z and i — z colours 
with axial ratio where the GZ spirals have been split into 
rough spiral types using the SDSS parameter /zjev (Figs 
B1-B3). This is described in Section 3.1.1 of the text; the 
g — z colour trends are shown in Figure 1101 



APPENDIX C: MONTE CARLO FITS: EXTRA 
FIGURES 

We show here in Figures [CHC2I the results for Monte Carlo 
realisations of the data split into subsamples by fDeV- De- 
tails are described in Section 4.1, and the results for the 
whole GZ spiral sample are shown in Figure [TS] 
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Figure Bl. The trend of u — z colour for GZ spirals separated into rough spiral types using the foeV parameter as indicated in the 
panels. As in Figure [9] the greyscale contours indicate the locus of galaxies, while overlayed are the median and interquartile range of 
data binned in log(a/b). 



© 2010 RAS, MNRAS 000, 1-25 



22 K.L. Masters et al. 



A 

i 1.0 



Pure disk, fdeV < 0.1 




0.0 



Disky, fdeV < 0.5 



Very bulgy, fdeV > 0.9 



0.0 



Bulgy, fdeV > 0.5 



0.0 0.2 0.4 0.6 
<- face-on log(a/b) 




0.0 0.2 0.4 0.6 0.8 1.0 
log(a/b) edge-on -> 



Figure B2. The trend of r — z colour for GZ spirals separated into rough spiral types using the fneV parameter as indicated in the 
panels. As in Figure [9] the greyscale contours indicate the locus of galaxies, while overlayed are the median and interquartile range of 
data binned in log(a/b). 
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Figure B3. The trend of i — z colour for GZ spirals separated into rough spiral types using the foeV parameter as indicated in the 
panels. As in Figure [9] the greyscale contours indicate the locus of galaxies, while overlayed are the median and interquartile range of 
data binned in log(a/b). 
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Figure CI. As Figure IT51 but only for galaxies with fneV < 0-1 ~ or "pure disk" spirals. Model parameters and Monte Carlo inputs are 
identical to that used for the full sample of GZ spirals. 
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Figure C2. As Figure IT51 but only for galaxies with foeV > 0.9 - or "very bulgy" spirals. Model parameters and Monte Carlo inputs 
are identical to that used for the full sample of GZ spirals. 
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